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In most introductory texts to quantum computation, such as |ChNi| 
or |Presj . the presentation of the quantum Fourier transform relies on the 
derivation of an explicit formula, which is then translated into a quantum 
circuit. The connection with the classical fast Fourier transform is also cited, 
but not explicited, or left as an exercise. 

In this note we construct a quantum Fourier transform circuit in a recur- 
sive way, by directly copying the construction of the fast Fourier transform 
algorithm that is given, for example, in |CLRSj . chapter 30. We do not 
pretend this presentation to be original, nor claim for any anteriority. The 
aim of this paper is purely pedagogical. 

The author would like to thank Prof. Hugues Randriam for his supervi- 
sion of the writing of this text. 

1 The fast Fourier transform 

Let n be an integer and N = 2 n . The discrete Fourier transform asso- 
ciates to each complex iV-tuple a = (ao, . . . , ojv-i) the complex iV-tuple 
b = (bo, . . .,b N -i) with 

JV-l 

(1) b k = 2"™/ 2 Cl k a 3 

3=0 

where 

(2) Cn = 
is a primitive iV-th root of 1. 
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Thus, apart from some power of v2, the discrete Fourier transform 
amounts to evaluating the polynomial 



N-l 

(3) P(X) = £ a,-*', 

3=0 

of degree 2 n — 1, at the powers of Q n . The fast Fourier transform algorithm 
performs this using a divide and conquer strategy. Obviously P can be 
written in a unique way 

(4) P(X) = P even (X 2 ) + XP odd {X 2 ). 

In this decomposition the degree of P even (resp. P oc id) is 2™" 1 — 1, and its 
coefficients are the dj for even (resp. odd) j's. Using the fact that = Cn-i, 
one gets 

(5) P(C) = WCti) + C^iWCti)- 

Thus, to perform a Fourier transform of order n, one has to perform two 
Fourier transforms of order n — 1, and then O(n) additions and multiplica- 
tions. We will summarize this as 

(6) FFT n = 2.FFT n _i + O(n) 

which finally gives us that the total number of elementary operations re- 
quired grows as 

(7) FFT n = 0{n2 n ). 



2 The quantum Fourier transform 

As we just saw, the exponential factor 2 n in the complexity of the fast Fourier 
transform algorithm comes from the factor 2 in equation ©, that is, from 
the fact that the computation of a Fourier transform of order n requires two 
Fourier transforms of order n — 1. Hopefully this factor 2 could be shrunk 
down to 1 using quantum parallelism. Indeed, recall that when one has a 
quantum circuit F acting on n — 1 qubits, then the circuit F Cg> 1 (see figure 
n below) on n qubits, applied to a vector \x) = |xo)|0) + gives as 

output (F|xo))|0) + (i ? |xi))|l), thus, in someway, computing simultaneously 
F\xq) and F\xi) at the cost of only one application of F. Carrying this out 
succesfully will lead to a dramatic decrease of complexity in the adaptation 
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n-1 




Figure 1: constructing F <8> 1 out of -F 



of the Fourier transform from the classical world to the quantum world: 
from an exponential complexity we will get to a quadratic one. 

To make it formally, consider an n qubits system, and denote the vectors 
of the computational basis by 

(8) \j) = \jl)-.-\jn-l)\jn), 

for j £ {0, . . . , N — 1} (recall that N = 2 n ) with base 2 expansion 

(9) J=jl2 n - 1 + ---+jn-l2 + jn- 

We will construct a quantum circuit which sends the unitary vector ^2,- a,j\j) 
to J2k bk\k) where the b^s are deduced from the a/s by relation (JJ). Remark 
that, since quantum circuits can only perform unitary transformations, this 
implies that formula Q is norm preserving. Although this fact is well known 
(Parseval's theorem), we will not need to assume it. Indeed our construction 
can be considered as a way of re-proving it. 
Now let 

N-1 

(10) |P> = J>,-|i>. 

3=0 

Grouping together terms according to j n = (that is, j even) or j n = 1 
(that is, j odd) we get a decomposition 

(11) |P) = |Pe«en)|0) + |P od d)|l), 

which is the exact analogue of decomposition 

According to formula (J5J), we have then to perform a Fourier transform 
on the first n — 1 qubits, which can be done as in figure ^ I n doing that 
one has however to be careful on the ordering of the output qubits. Indeed, 
write 

(12) k = ki2 n ~ l + k 2 2 n ~ 2 + ■■■ + k n ^2 + k n . 
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Then, as Cn-i = 1 and Cn™ = ~~ 1; formula |JSJ) can be re-writen as 

(13) P(C) = PeveniCtl) + (~ ^ C% Podd(d) 

where 

(14) A/ = k 2 2 n ~ 2 + ■■■ + k n ^2 + k n . 

Thus, while \P eve n) and \P dd) are vector states of the first n — 1 qubits 
ji, . . . , j n -i, their Fourier transforms should be output on the last n — 1 
qubits A>2, • • • , k n . To circumvent this difficulty we will consider a modified 
quantum Fourier transform, mQFT, which reverses the order of the output 
qubits. Obviously this is not a serious problem, since the correct order can 
then be restablished using only n/2 swap gates. 

Our modified quantum Fourier transform circuit of order n should then 
look like: 
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Figure 2: recursive scheme for mQFT, 



At point A we input the vector 

(15) \P} = \P even )\0) + \P odd }\l), 
so that at point B we get 

(16) |P^)|0} + |P^d)|l), 

where we denoted by 7 the effect of our (modified) quantum Fourier trans- 
form. According to formula (|13|) we have 



(17) P{£) 



Peven(^'-l) + ^ PoddiQ^l ) ^ h = 
Peven(C'-l) ~ C£*WC£-l) for h = 1 
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so that, putting the 2 n l 2 back, at point C we should get 



(18) \P) = ^={\Peven) + g\Podd))\0) + ^=(\P. 



Cn fc 'l^))|l) 



where Cn can be understood as the diagonal matrix acting on the first n — 1 
qubits by multiplying \ k') by , or equivalently as the controlled operation 
that multiplies the last qubit by ^ when the first n — 1 qubits are set to 
\k>). 

Clearly one can go from (fTHf) to (fTHj) bv sending |0) to ^(|0) + |1)) and 
|1) to -^(Cn |0) ~~ Cn |1))> that is by acting on the last qubit through 



(19) 



1 

72 



i 



it' 



where H = ^ 



U P^ 2 p^n— 1 pfen 

is the usual Hadamard matrix and i?^ 



2l7T 





1 1 
1 -1 

This means that we can complete our quantum Fourier circuit by replac- 
ing the '?' circuit in figure |2] with 



1 
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Thus we see that the (modified) quantum Fourier circuit of order n can 
be constructed from the order n— \ one by adding n — 1 controlled operations 
and one Hadamard gate. Putting all together, since the quantum Fourier 
circuit of order is trivial (the identity) , we get a total complexity of ra ^ n 2 t " 1 ^ 
gates. 

The reader will remark that the circuit we get is the inverse of the one 
that is given in gure 5.1. This is not a mistake: it only indicates 

that the Fourier transform is essentially its own inverse. 
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